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... A suggested solution to the dark energy problem is the void model, where accelerated expansion is 

CO ' replaced by Hubble-scale inhomogeneity. In these models, density perturbations grow on a radially 

T-H I inhomogeneous background. This large scale inhomogeneity distorts the spherical Baryon Acoustic 

, Oscillation feature into an ellipsoid which implies that the bump in the galaxy correlation function 

occurs at different scales in the radial and transverse correlation functions. We compute these for the 
5_( ' first time, under the approximation that curvature gradients do not couple the scalar modes to vector 

^ , and tensor modes. The radial and transverse correlation functions are very different from those of 

' the concordance model, even when the models have the same average BAO scale. This implies that 

' if void models are fine-tuned to satisfy average BAO data, there is enough extra information in the 

T-H , correlation functions to distinguish a void model from the concordance model. We expect these new 

OA ' features to remain when the full perturbation equations are solved, which means that the radial and 

transverse galaxy correlation functions can be used as a powerful test of the Copernican Principle. 

O: 

^ ; I. INTRODUCTION 

O 1^ The homogeneous and isotropic Friedmann-Lemaitre-Robcrtson- Walker (FLRW) universe is based on the assump- 
Q ' tions that the Copernican Principle holds true and that the universe is isotropic about our location. The high degree 
5^ ' of isotropy of the Cosmic Microwave Background (CMB) provides strong support for the second assumption. The 
. Copernican assumption cannot be directly tested, although null tests of this assumption have been devised (for re- 
d ' views, see [HQ)- A spatially flat FLRW universe containing 20% cold dark matter (CDM), 5% baryons and 75% dark 
energy in the form of the cosmological constant (A) provides an excellent fit to the wide range of observations to date. 
■ However, there is still no satisfactory theoretical explanation from fundamental physics for the observed value of A. 
J> ' This has prompted many authors to consider alternatives to ACDM, such as modified gravity Q, backreaction and 
^sj • spherically symmetric but inhomogeneous exact solutions to Einstein's equations known as Lemaitre-Tolman-Bondi 
O; (LTB) models (see 0,3 for recent reviews). 

^ : In this work we focus on the latter alternative. Hubble-sized LTB void models break the Copernican Principle 
' I ' by placing our galaxy in a special position, at the centre of an underdense region of O(Gpc) scale. The simplest of 
\^ . these models are able to fit the distance-redshift data from type la supernovae (SNIa) and the CMB, though this 
' gives considerable tension with the locally measured value of Hq (see e.g., @). Attempts to overcome this require the 
[ relaxation of assumptions on the homogeneity of the early universe in the form of isocurvature modes or a change 
to the primordial power spectrum Q. 

It remains an open issue whether they can fit the Baryon Acoustic Oscillation (BAO) data, and the more general 
data on the growth and distribution of large scale structure. Previous papers [ol-ll^ have computed the BAO scales 
in a geometric approximation, using the anisotropic expansion rates of the background model but ignoring any effects 
, from the anisotropic growth of structure in void models. This is not surprising, since structure formation on an LTB 
background has not yet been calculated, because it is much more complicated than in the standard model p^l - ITsj . 

Here we calculate for the first time the 2-point correlation function on an LTB background, and use it to extract 
the radial and transverse BAO scales. This incorporates the effects of the evolution of density perturbations on an 
LTB background, using the perturbation formalism developed in [l8| . We neglect the coupling of scalar to vector and 
tensor modes in the metric potentials. This is expected to be a good approximation for the simplest LTB models, in 
which the background shear — responsible for the coupling of the modes in the first place — is typically of the order 
of a few percent (The accuracy of this approximation is under investigation via numerical solutions [Tgj.) 

LTB models have enough freedom to always fit the average BAO scale. In these models the proper radius of the 
sound horizon at the drag epoch is approximately given by 0, 0] 
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where A^eff ~ 3.04, ff, ~ flf,/^m is the local baryon fraction, rj = 10 ^'^tjiq is the baryon-photon ratio at that time. 
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Td is the temperature at the drag epoch and it is assumed that during the process of recombination, the scale of 
the void inhomogeneity is much larger than the horizon size at that time (~ 100 Mpc). In general, /f, and r^io have 
radial degrees of freedom in them, and are no longer measured by the CMB at the radial scales of interest for the 
BAO. Consequently, ds is not constant spatially, and can vary over the scale of the model. Therefore, measuring the 
mean BAO scale in some shell around us cannot place constraints on late-time inhomogeneity without some other 
measurement of fb and r] in the same shell at early times - which lie inside our past lightcone. In fact, the models can 
even be fine-tuned to have the same radial and angular BAO scales, by altering the bang time function appropriately. 

As we shall show, however, the radial and transverse 2-point correlation functions contain much more information 
than the peak positions which determine the BAO scales. In real space, the radial and transverse correlation functions 
are typically very different from each other. To make them close to those of the ACDM model with the same BAO 
scale would require high levels of fine-tuning in either dg or the primordial power spectrum. Effectively, the radial 
and transverse 2-point correlation functions can thus be used as an important probe of the Copernican assumption. 

This paper is organized as follows. In Section |lT] we recap the standard geometric approximation for computing 
BAO scales and describe our method for determining the background dynamics of the LTB spacetime. In Section 
mil we provide an overview of perturbation theory in LTB via the 2+2 decomposition approach. This is followed by 
a derivation of the anisotropic two-point correlation function for the gauge-invariant matter density perturbation in 
Section HVl where we present the computation of the correlations and the BAO scales. Finally, in Scction|V]we discuss 
the consequences of our results. 



II. BACKGROUND MODEL AND EVOLUTION OF THE BAO SCALES 



The background void model is described by the LTB metric, 

a\it r) 

ds^ = -dt^ H L_— (ir^ + a\ (t, r)r^dn'^ , a» = ^a' , ^a raj_ , (2) 

1 — K{r)r^ " 

where c?a is the angular diameter distance, and a prime indicates d j dr. In the FLRW limit, a|| (i, r) = aj.(t, r) = a(t) 
and nir) ~ K . The expansion rates transverse to and along the line-of-sight are 

H^{t, r) = — — — and H\\ (t, r) =^ " . (3) 
a^(t,r) o,\\\t,n 

The past lightcone of the central observer has null geodesies that are given by 

dt _ 1 dr _ [1-K(r(z))r2(z)]'/' 



dz (l + z)i?||(t(^),r(z)) ' dz {l + z)a\\{t{z),r{z))H\\{t{z),r{z))' 



(4) 



We use the notation F{z) = F{t{z), r{z)) to denote evaluation on the past lightcone. 

The anisotropic expansion rates ([3]) act on the acoustic sphere of proper radius L* at an initial high redshift z*, so 
that by redshift z it has evolved into an axisymmetric ellipsoid with semi-axes 



II a||(t*,r(z)) a±(tt,r{z)) 

However, this geometric approximation does not give the correct BAO feature in the galaxy distribution ~ because 
it neglects the anisotropic effects of perturbations in LTB and their impact on the correlation function. Previous 
work 0, H, [lol - [l2i [2l| on comparing the BAO scales in LTB with observations has all neglected the effects of LTB 
perturbations. Below we fill this gap by computing the correlation functions associated with the density perturbation 
and then extracting the BAO scales from the correlation functions. 

The observable quantities of the BAO feature are its redshift extent 5z{z) and angular size S9{z). These are 
converted to the physical radial and transverse length scales via 

H^)= V L^=dA{z)5eiz), (6) 

{l + z)Hn{z) 



for small 5z and 66. Note that we neglect redshift space distortions. 
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The quantity Q 
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encodes an average of the two observable scales of the sound ellipsoid. In an FLRW model it reduces to 

1/3 



L*(l + 2*) , . 

dz = , Dv{z) = 
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where L^,{1 + z*) is the comoving sound horizon and Dy is the standard volume-averaged BAO scale [22| . 
The LTB analogue of the Friedmann equation is 



Hioir) aiit,r) ai(t,r) 



, where n„,{r) + ^Kir) ^ 1, a^{ta,r) = 1, H^r) = H^{to,r). (9) 



The observed Hubble constant is iJp" = H±o{0) = lOO/i km/s/Mpc, where "in" indicates evaluation at the centre. 
For open LTB models, the parametric solution is 



a±{t,r) 



il,n{r)[cosh2u{t,r) — 1] 

2[l-n^{r)] ' 
nm{r)[sinh2u{t, r) -2u{t,r)] 



2H^o{r) [1 - nm{r)f^ 

where we choose a simultaneous big bang (uniform bang time function). Setting t ~ in (jlO[) and (jlip gives 



- 1 



r2,„(r) [sinh2ito(r) — 2ito(r)] 



2tQ [1 -f7„(r)] 
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(10) 
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Thus H±o{r) is determined when D,m{r) and to ^-rc chosen. Then (jlip determines u{t,r) and a^(t, r) follows from 

m- 

For the purposes of this study, wc choose a simple Gaussian void profile for the dimensionless density parameter, 



(13) 



n^ir) = nT - i^T ~ f^™) exp , with = i, = 0.7 , 



where "out" refers to the asymptotic Einstein-de Sitter region, and a characterizes the size of the void (see [2J] for 
details). The physical matter density is then 



8TrGan{t,r)a\{t,r) 



3 + r <^ 2 



(14) 



III. SCALAR PERTURBATIONS ON AN LTB BACKGROUND 

The full perturbation theory on an LTB background is developed in [l^ via a 2+2 split of the spacetime, which 
makes explicit the coupling of vector and tensor modes to scalar modes at linear order. A first approximation is to 
neglect this mode-mixing, and focus only on 'scalar' modes which occur in the even parity sector. Then the perturbed 
metric in Reggc- Wheeler gauge is ([l^: with notational change, Lp — > —2$) 

ds^ = - [1 + 2$(i, a;)] dt^ + [1 - 2$(t, x)] g.,jdx'dx' , (15) 

where gij is the spatial part of The Newtonian potential obeys a simple generalization, without gradients, of the 
FLRW evolution equation for the Newtonian potential [Tsj : 



$ + 4ff^$ - = 
ai 



(16) 
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FIG. 1: Upper: Background density and expansion rates for the models H49 [) - (|5H) . Lower: Using the geometric approximation 
Q , the evolution of the BAO length scales (left), and the average BAO scale (right). Black circles indicate measurements from 



Because there are no spatial gradients, $ evolves independently in each r = const shell, as if in a separate dust 
FLRW model. This does not mean that there is no dependence on spatial gradients: density fluctuations depend on 
spatial gradients of $ which couple to the anisotropic expansion of the model. The gauge-invariant matter density 
perturbation A is found via the equivalent of the Poisson equation in LTB : 



AirGajpA = £ [$] , 
where £ = (1 - Kr^)d^ 
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a|(F|| + 2Hi_)dt - a\H^{Hi_ + 2H» 



(17) 



(18) 



(Recall that the LTB model contains only CDM and baryons.) In FLRW, we recover the standard Poisson equation: 



3K 



$ - 3a2i7($ + H$), 



where = (1 - Kr^)d^ 



(2 - 3Kr 



-dr — 



£{£ + 1) 



Here £ is the angular wave number in a spherical harmonic expansion, 

$(t,a;) = ^$£,„(t,r)r^™(6l,^). 



(19) 
(20) 

(21) 



and similarly for A. 

We set initial conditions for $ at a high redshift, = 100, where we assume the background is effectively FLRW. 
We write 



$£m(i,?-) = <^(t,r)$,f™(r), (l){U,r) = 1. 



(22) 
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FIG. 2: The gravitational potential (/> as a function of radius today (left), and of redshift (right). 
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FIG. 3; The normalized density perturbation [(i?o'")^Afm/"I>t^m| against redshift, for small (left) and large (right) i. 



The subsequent time evolution of (/)(i, r) is then determined by (|16p for each r. Using (fTO]) and (jlip . this implies 

(j){t,r) = C{r)- — g-^^[sinh2ii(t,r) - 6u(t,r) +4tanhu(t,r)] , (23) 
sinh u{t, r) 

cosh u* (r) [sinh 2u*(r) — 6w*(r) + 4 tanh (r)] 



Now, ^*irnir) can be written as 



^*im{r) = \l-i' I d^kMkr)<^,{k)Yemik), (25) 



which is related to the power spectrum via 

27r2 



($* (fei) (fc2)) - -rr^** (fci) ^ (fci + '^2) . (26) 
The initial power spectrum of the Newtonian potential is given by 

r^,{k)^^rn{ko)T^{k), (27) 
25 

where Vn{ko) = 2.41 x 10~^ is the amplitude of the primordial curvature perturbation on the scale k^ = 0.002 Mpc~^, 
and T{k) is the matter transfer function, with T(ko) ~ 1. The concordance parameters ()49p are used in the fitting 
formula of [25| to compute r(fc), which is employed in all of the models. 
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Note that when using a flat FLRW initial power spectrum, wc need to use the flat FLRW comoving coordinate rp 
in (|25p at t^, as opposed to the LTB coordinate r. Proper radial distance is independent of coordinates: dp(t^,,rp) = 
dp(t^,,r). Since dp{t^,,rp) — a(t^,)rp, we find that 



rp = {l + z.) I dr^h^ ^ f{r), 
^l-K(r)H 



where /(r) w (1 + z»)aj^(t», r)r since -^/l — K{r)r'^ ?» 1 for all r and a|| = dr{a±r). Then becomes 

- f d^kMkf{r))<^>4k)Ye^{k). 



(28) 



(29) 



IV. CORRELATION FUNCTIONS AND THE BAO SCALES 

The two-point correlation function (2PCF) for the density perturbation A, as observed by a central observer down 
their past lightconc, is 



C(ii, ^2, r-i, ra) = (A(<i, ri)A(t2, rz)) = ^(^(^i), ^(^2), rizi)fi,r{z2)r2) 
= ^{t{zi),t{z2),r{zi),r{z2),S9), where ri • r2 = cos (50. 



(30) 



The second line follows from statistical isotropy, which applies for central observers. We neglect redshift space 
distortions for simplicity, since wc are not testing the void models against data but only comparing them with the 
concordance model. We also neglect all complications from bias for the same reason. 
Using the Poisson equation ([T7)) . the correlation function (PH]) becomes 

£,{zi,Z2,6e) = [(47rGa||^a||2)^PlP2] ^ ^ /^l^l 'C2?i2>£m(ri)Y«'m'(r2)($*l£m$*2f'm') 

= [87r3(Ga||^a||2)ViP2]"' ^ ^ d^fci d'fe A [0ij£(fci/i)] £2 [02j£' (^2/2)] >£m(ri)y£'™'(f2) 



X ($*(fcl)*,(fc2))y£m(fcl)yrm'(fc2), 



(31) 



where a subscript i = 1,2 on a function of (t,r) means the quantity is evaluated at {t{zi)^r{zi)). Using (|26p and 
standard identities in (|3T|) . we get 



/J 1 
-^Mzi,k)JpXz2,k)V^{k), (32) 

(33) 



where J«(z,fc) = C[4>{t{z),r{z))j,{kf{z))\ 
To evaluate (|33)) . we use ((T8)) and the following identities for the spherical Bcssel function 



The result is 



where 
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(40) 
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and 

(42) 

u,_L L \ / 

In the flat FLRW case, ([M]) becomes 

j£(z,fc) = - ■ia^H^+{2,c?B'^ ^k^^<i) ji{kr) . (43) 
The usual flat FLRW correlation function may then be obtained using the following identity: 

^(2£+ l)P,(cos(5^^)J,(fcrl)J,(fcr2) - , (44) 

f=0 



where s = s{zi,Z2,S9) = \/r\^r\~ Irxr^ cos bd. 

In LTB, the real-space radial and transverse BAO scales are different, and are given by the peaks in the radial and 
transverse correlation functions. These we define as: 



radial2PCF: (z, Jz) = ^(^1, ^2, 0) = ^(2^ + l)c|' (z, (5z), (56* = 0, z = zi, = za - , (45) 



transverse 2PCF: C^(z,(56') = i{z,z,b6) = ^(2^ + l)PKcos 56l)C/(z), (46) 

I 

where the radial and transverse coefficients follow from (|32]) : 



1 C r\h 

(z,<5z) = [(4^Ga||^a||2)ViP2] — Ji{z,k)Je{z + Sz,k)V^{k), (47) 
C,^(z) = (AnGajpY' 1 J,\z,k)r^{k). (48) 



Equations (|32l) . ()36p and (|45 |) -(j48 |) summarize our new results that derive the correlation function of matter density 
perturbations on a radially inhomogeneous background. 

A. Computation of the correlation functions and extraction of the BAO scales 

We can now compute the correlation functions for two specific LTB models, and compare with the standard case 
(see e.g. [26l - l3l| for various approaches to compute these quantities from galaxy surveys in the standard homogeneous 
framework.) We consider the following models: 

• Flat ACDM: a concordance model, with ^bh^ = 0.02273 and VLch?' = 0.1099, as given by WMAP 5-year 
CMB-only best-fit results (see Table 6 of [11]). Setting h = 0.7, this implies 

fACDM: VL,n = VLb + VLc^ 0.2707, f^A = 1 - = 0.7293 . (49) 

This is our benchmark model which we use to compare with void models of the type given by (|13p . 

• SV: a small void (compared to those that fit SNIa luminosity distances [l^]) of type ([13]), with 

SV: VlZ = Q.2, CT = 500Mpc. (50) 

• BV: a big void of type p^ . chosen so that its anisotropic expansion rates provide a good fit to observations of 
the average BAO scale (O. We performed a fit to measurements of (see Table 3 of [l^l), and found the 
following best-fit parameters: 

BV: rjj^J = 0.32, o- = 4.84 Gpc. (51) 
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FIG. 4: Angular power spectra + 1)C/ at various redshifts for the BV and fACDM models. (The drop in power on small 
scales is due to smoothing the power spectrum below 1 Mpc.) 



,x 10 



3 Radial 2PCF:z^ =0.05 




0.02 0.04 0.06 0.08 
8z 



,x 10 



Jransverse 2PCF: Z|=0.05 




20 40 60 

86 (degrees) 



FIG. 5: Radial (left) and transverse (right) correlation functions for the SV and fACDM models at zi — 0.05. 



In the two void models, we choose FLRW initial conditions to ensure that the effects we find arise from the evolution 
of structure on the inhomogeneous background. We take the early-time parameters fb and in ([1} to be those derived 
from the same WMAP 5-ycar values used for the fACDM model. This fixes the initial proper BAO scale to be the 
same in all models. The background density flm and expansion rates , H± are shown for these 3 models in Fig. [T] 
(upper panels). We also show (lower panels) the geometric approximations to the radial and transverse scales, -^l''" 

and L^°, and the average BAO scale calculated from them. 
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Figure [2] shows the current profile and the redshift evolution of the gravitational potential for the 3 models. Note 
the greater decay in the amplitude of 4> for the void models, due to the presence of curvature, which explains the 
decrease in the overall amount of clustering relative to ACDM. The normalized density perturbation is illustrated 
in Fig. [3] for the 3 models. For small-scale modes (large £), A scales approximately as (1 + z)~^. For the large-scale 
mode £ = 2, the 'decaying' behaviour at high redshift is due to the mode entering the Hubble-scale at low redshift. 




FIG. 6: Radial (upper) and transverse (lower) correlation functions at various redshifts for the BV and fACDM models. 

We calculate the correlation functions by smoothing away power on scales below IMpc, via 7^$(fc) 
7-'$(fc) exp[— fc^/(l Mpc~^)^]. This makes the sums over £ in the correlation functions ps)) . (|46p converge relatively 
quickly (typically we require £nia,x{z) ^ 10r(z)/Mpc), but without altering the resulting correlation function. Figure|3] 
shows the angular power spectrum £{£ + 1)C^ for the BV void model compared to the concordance model. The drop 
in power for high £ results from the smoothing. 

The correlation functions for the two void models arc shown in Figs. [5] and [HI The radial correlation 

function ^H, starting at various redshifts zi and extending to Z2 = zi + Sz, shows the correlation of structure along a 
line of sight, as the observer looks into higher density regions. The redshift extent of the radial BAO feature is (Szpoak, 
which is given by the location of the bump in ^'L The transverse correlation function describes the correlation 
across the sky in a sphere at redshift zi. The angular size of the BAO is (50pcak, given by the bump in 

It is apparent from Fig.[5]that for the SV model, the radial correlation function is very different from the concordance 
one. This is due to the large curvature gradients at low redshift, compared to void models that fit SNla data. We 
neglect redshift space distortions - but curiously, the effect of the void is qualitatively similar to the effect of redshift 
space distortions in FLRW (see [33j). 

We determined ^Zpcak and S9peak numerically from the local maxima in the correlation functions. The results are 
shown in Figs. [7] and |S1 In these figures we also show the geometric approximations (i.e., without incorporating the 
effect of perturbations) , 

,5zs™ = Lf°(l + z)Hn, 69^°° = — i-, (52) 
" dA 

where L^\''°, L^±° are given by ([5|). Our results show that the geometric formulas commonly used for constraining LTB 
with BAO fail at the percent level. While current data are not able to resolve such differences, this may be possible 
with future surveys such as SKA and Euclid. Furthermore, note that the size of these corrections are of a similar 
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order to the corrections from redshift space distortions in FLRW [33]. Note that the geometric formulas in (|52p give 
the correct observed scales for fACDM - except for large 50, for which the small-angle formula in ([5^ breaks down. 




FIG. 7: Upper: Redshift extent (left) and angular size (right) of the BAO feature for the SV and fACDM models. Lower: Ratio 
of the BAO length scales. In all plots, we show the results of the full calculation based on the correlation functions (numerical) 
and of the simplified geometric approximation. 



V. CONCLUSIONS 



We have derived for the first time the anisotropic real-space two-point correlation function for the gauge-invariant 
matter density perturbation, in an LTB universe with radial inhomogencity in the background - summarized in (|32l) , 
([55)1 and (|l5| (|15|) . For this we neglected the coupling of scalar modes with vector and tensor modes - which should 
be a good approximation, at least on the large scales relevant for the BAO. An analysis of the effects of mode-coupling, 
which would entail the integration of partial differential equations, is currently underway [l^. We also neglected bias 
and redshift space distortions, since our primary focus was a comparison with the concordance model, not to test 
void models against data. Redshift space distortions in LTB void models deserve further investigation, in particular 
to check whether the FLRW formula provides a useful approximation. 

We computed the radial and angular correlation functions for two void models, one relatively small (SV) and one 
Hubblc-sizcd void that fits the average BAO data (BV) - see Figs. [S] and [S] We used the peaks of the computed 
correlation functions to extract the radial and transverse BAO scales. The results were compared with the geometric 
approximation that has been used in all previous work, showing that the geometric approximation to the BAO scales 
in LTB fails at the percent level - see Figs. [7] and HI Future large- volume surveys, such as SKA and Euclid, may thus 
be able to rule out the void models on the basis of their BAO scales. 

However, even if void models can be fine-tuned to reproduce the radial and transverse BAO scales, these scales 
represent only one feature in the galaxy correlation functions. The void correlation functions differ significantly from 
those of the concordance model (Figs. [5] and [5]). In particular, the void radial correlation can become negative (anti- 
correlation) before, and even at, the BAO peak, where the concordance correlation is positive. The void transverse 
correlation may be positive for all scales, unlike the concordance one. These features resemble the effect of redshift 
space distortions in FLRW (see Figs. 4 and 6 in [s^), since the anisotropic expansion rate in LTB can mimic the 
effect of radial peculiar velocities in FLRW. However, there are significant further differences between the two models 
which arise from the effect of LTB perturbations. 
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FIG. 8: As in Fig. [3 for the BV model. 



This leads to our key final result: even if the radial and transverse BAO scales match observations, the radial and 
transverse correlation Junctions contain direct signatures of the anisotropic growth of perturbations in a non-FLRW 
model. These correlation functions can thus be used as direct tests of the Copernican Principle. 
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